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The transmission properties of armchair graphene nanoribbon junctions between graphene elec- 
trodes are investigated by means of first-principles quantum transport calculations. First the depen- 
dence of the transmission function on the size of the nanoribbon has been studied. Two regimes are 
highlighted: for small applied bias transport takes place via tunneling and the length of the ribbon 
is the key parameter that determines the junction conductance; at higher applied bias resonant 
transport through HOMO and LUMO starts to play a more determinant role, and the transport 
properties depend on the details of the geometry (width and length) of the carbon nanoribbon. In 
the case of the thinnest ribbon it has been verified that a tilted geometry of the central phenyl 
ring is the most stable configuration. As a consequence of this rotation the conductance decreases 
due to the misalignment of the pi orbitals between the phenyl ring and the remaining part of the 
junction. All the computed transmission functions have shown a negligible dependence on different 
saturations and reconstructions of the edges of the graphene leads, suggesting a general validity of 
the reported results. 
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I. INTRODUCTION 



Graphene, a perfect carbon monolayer 

sp^— hybridized, has attracted a huge interest since 
its discoveryJ^ Besides a pure theoretical interest, 
its possible applications in carbon-based electron- 
ics represent a very exciting perspective. Graphene 
displays a very peculiar electronic structure, arising 
from the confinement of electrons in two dimensions 
and its geometrical symmetries. Graphene is a zero 
gap semimetal whose specific linear electronic band 
dispersion near the Brillouin zone corners (Dirac points) 
gives rise to electrons and holes that propagate as 
massless fcrmionsj^"— Graphene nanoribbons (GNRs) 
are one dimensional graphene strips that are considered 
as promising candidates building blocks for future 
electronic applicationsi^"— Several methods allow the 
production of GNR, including mechanical cutting of 
exfoliated graphene;^ patterning of epitaxially grown 
graphene^ bottom-up chemical mcthodsfii and carbon 
nanotubes unzippingJ^ Due to the honeycomb geometry 
of graphene, GNR can be patterned along two prefer- 
ential directions giving rise to armchair shaped edges 
graphene nanoribbon (aGNR) and zigzag shaped edges 
one (zGNR). The finite size of the GNR gives rise to 
a large variety of electronic behaviors that could be 
relevant in transport. Considering the case of aGNR, 
it has been demonstrated that its electronic properties 
sensibly vary by changing the width of the ribbonJ^^— 



In fact, the width and the electron energy gap A are 
related to each other primarily in inverse proportion. In 
particular, there exist three different classes of iV-aGNR 
(aGNR with N dimer lines) for what concerns the value 
of the gap Aat: Asp+i > Asp > Asp.i, with p integerii^ 
If we restrict to each of these classes, the energy gap 
decreases as p increases. Transport properties of pristine 
GNR have been studied both for the ideali2r— and 
defective^Sr— case. Electronic transport has also been 
studied, by a tight-binding approach, for junctions 
connecting zGNR of different widths^^"— revealing 
the crucial role played by corner edge structures)^ A 
challenging technological issue would be to exploit the 
high mobility properties of graphene and the finite-size 
characteristics of GNR, by creating heterostructures 
which exploit both these advantages. These systems may 
be used as interconnections to transmit signals in future 
pure-C-based electronic devices. In fact, the efi^ects of 
quantum confinement in graphene nanoconstrictions 
have been studied showing their analogy with optics 
phenomena^ as well as their exploitation for valley filter 
applications;^ and recently a quantized ballistic con- 
ductance has been measured in such structures »22i Even 
more complex junctions have been experimentally built, 
with carbon nanotubes as interconnections between 
graphene bilayersj^S However, due to their fiat structure, 
GNRs seem easier to pattern than carbon nanotubes. 
Owing to the large zoology of hybrid graphene inter- 
connections, it is of great importance to perform a 



systematic investigation of their transport properties at 
the simplest level of configurational complexity, in order 
to clarify how the basic geometrical parameters affect the 
conductance. Most of the previous works investigating 
GNR junctions and graphene nanoconstrictions rely on 
standard tight-binding calculations. It has been shown 
by Louie and co- workers^ that the predictions of simple 
tight-binding models on GNR may lead to incorrect 
band structures and energy gaps, since the bonding 
characteristics between atoms substantially change at 
the edges. An ab-initio study of GNR junctions avoids 
such these difficulties, providing an accurate description 
of their transport properties^, which are directly related 
to the underlying electronic band structure. 

In this paper, we present a first-principles study, by 
means of the non-equilibrium Green's functions (NEGF) 
technique, of the electronic and transport properties 
of systems consisting of two semi-infinite graphene lay- 
ers interconnected by an hydrogen-passivated armchair 
graphene nanoribbon. The transport properties of such 
junctions are predicted to strongly depend on the GNR 
geometry, while are quite robust to changes of the 
graphene's edges geometry. These structures combine the 
high mobility of graphene electrodes^"— with the intrin- 
sic semiconducting behavior of GNR. We show that these 
semiconducting hybrid graphene-GNR junctions have a 
significant gap in the transmission spectrum, which may 
be exploited to build logic devices which require a large 
on-ofF ratio in the current i^ In the following discussion, 
we will show the properties of aGNR junctions consid- 
ering one representative for each class according to their 
width: we will take into account 3-aGNR, 5-aGNR, and 
7-aGNR. In the case of 3-aGNR, we will consider rib- 
bons of three different lengths, namely those consisting 
of four (4L), six (6L), and eight (8L) zigzag lines along 
their axis. The 8L 3-aGNR will be considered as an il- 
lustrative model to study the effects of the application 
of a bias. We will discuss the effect of the rotation of a 
phenyl ring, and show how different edges configurations 
affect the transport properties. 




FIG. 1. System setup for 4L 7-aGNR. The left and right 
electrodes and the extended-molecule region are highlighted. 



plemented in SIESTA packagei^i^ The exchange- 
correlation energy and electron— ion interaction are de- 
scribed by the Perdcw-Burke-Ernzerhof (PBE)^ gen- 
eralized gradient approximation (GGA) and norm- 
conserving pseudopotentials^ in the fully nonlocal form, 
respectively. A double-C polarized basis set of numeri- 
cal atomic orbitals is used and the energy cutoff for real- 
space mesh is set to 200 Ryi^ Preliminary tests indicated 
that the relaxation of the carbon atoms in the leads did 
not affect the transport properties of the systems under 
study, so in most cases we considered non-relaxed geome- 
tries. The edges of semi-infinite leads are saturated with 
one relaxed hydrogen per carbon atom. We also verified 
that the relaxation of the aGNR docs not affect signifi- 
cantly the electronic and transport properties of the sys- 
tem. For the calculation of the transmission coefficients, 
60 fc-points along the transverse direction in the 2D first 
Brillouin zone are used. Periodic images of the graphene 
layer are separated by 15 A along the normal direction. 
The electronic transport is studied with the TranSIESTA 
code^^ which combines the NEGF technique with den- 
sity functional theory. The transmission function of the 
system can be obtained by the following equation: 



II. MODEL AND METHODS 

The system setup for a 4L 7-aGNR junction is shown 
in [T] as a representative case. This open system is con- 
stituted by three parts: the left (L) and right (R) semi- 
infinite graphene leads, and an extended-molecule (EM) 
region. The junction is constructed so that the periodic 
replicas of the aGNR along the direction parallel to the 
electrode edge are separated by 7.43 A. This corresponds 
to three unit cells of graphene and we verified that the 
interactions among the replica are negligible. The width 
of the leads part included inside the EM region is chosen 
after converging the transmission function: six carbon 
zigzag lines per side are sufficient to reach stable results. 

The electronic structure calculations are carried out 
using the first-principles self-consistent method im- 



T{E, V) = Tt[Ti^{E, V)G{E, V)Tn{E, V)G^{E, V)], 

(1) 
where the spectral density Fljr) describes the coupling 
between the L (R) electrode and the EM region and it is 
given by the imaginary part of the electrode self-energy: 
rL(R)(^,^) = *(Sl(R) - S*(j^))/2. The self-energy de- 
scribes the hopping across the surface separating one 
lead and the EM region, and establishes the appropriate 
boundary conditions for the Green's function calculation. 
G is the retarded Green's function of the EM region, for- 
mally given by: 



G{E, V) = [ES - H{V) - Sl(£;, V) - ^^{E, V)] ' , 

(2) 
where 5* is the overlap matrix and H is the Hamiltonian 
of the system when a bias voltage V is applied. The 




FIG. 2. Electronic and transport properties of 4L 3-aGNR. 
Upper panel (a): DOS (solid line) of the junction and PDOS 
(shaded grey area) on the ribbon region. Red triangles repre- 
sent the eigenstates of the isolated (unconnected and hydro- 
genated) 4L 3-aGNR linker. Lower panel (b): transmission 
function. 



current is simply given by the following integral: 

I{V) = — y AE T{E, V) [f{E ~ a^l) - !{E - m)\ , 

(3) 
with / being the occupation Fermi function, /iL(R) the 
chemical potential of the L (R) electrode, and V = [-L\,— 



III. RESULTS AND DISCUSSION 

We first investigate in detail the electronic properties 
and the transmission function of the 4L 3-aGNR junc- 
tion. For the other systems, similar considerations can 
be done. In[2]we show the density of states (DOS) in the 
EM region and the projected density of states (PDOS) 
on the ribbon region for the 4L 3-aGNR, in comparison 
with the transmission function. Around the Fermi en- 
ergy {Ep), the DOS resembles that of graphene, in fact 
it goes to zero almost linearly with a deviation due to 
the presence of vacuum portions in the molecular bridge 
region. The signature of the molecular energy levels is 
clear both in the DOS and PDOS, and the transmis- 
sion function has higher intensity in correspondence to 
those peaks. We also consider the isolated molecule ob- 
tained by cutting the bonds between the 4L 3-aGNR and 
the graphene leads, and saturating them with hydrogen 
atoms. The red triangles represent the eigenvalues of the 
isolated molecule, which correlate well with the position 
of the resonances in the DOS, PDOS, and transmission. 
The (P)DOS for the other studied systems have a similar 
behavior, and are not reported here. 
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FIG. 3. 4L 3-aGNR junction: isosurface plot of the wave- 
functions at energies corresponding to the LUMO and HOMO 
resonant transmission peaks, calculated at fcy — 0. 
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FIG. 4. Transmission function of 3-aGNR junctions for three 
different lengths: 4L (red), 6L (green), and 8L (blue). The 
dotted line is the transmission function of the pristine infinite 
3-aGNR. 



Moving away from Ep, the transmission function in- 
creses up to significant values in correspondence to the 
two peaks at —1.2 and 1 eV; these peaks are generated 
by the hybridization of the molecule's highest occupied 
molecular orbital (HOMO) and the lowest unoccupied 
molecular orbital (LUMO) states with the leads, play- 
ing the role of channels between the two graphene elec- 
trodes. The corresponding wavefunctions of the interact- 
ing system (shown in [3]) calculated at /c|| = 0, have the 
same shape and symmetry of that of the isolated molecule 
(see [8]). These wavefunctions represent good conducting 
states, as they are dclocalized along the molecule and 
propagate inside the electrodes with the same symme- 
try. Some other states farther from Ep contribute to the 
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FIG. 5. Transmission function of 4L A''-aGNR junctions for 
A^ = 3, 5, and 7. Dotted line: transmission function of the 
infinite pristine A^-aGNR. 



transmission. All these states are generated from the tt 
orbitals of the C atoms and they are delocalized along 
the junction. 

We now discuss the dependence with the length of the 
ribbon forming the junction. The transmission functions 
for the 4L, 6L, and 8L 3-aGNR junctions, reported in 
m show some general trends. The peaks become denser 
when the length of the junction increases, because they 
are related to the discrete structure of the electronic 
states of the nanostructured ribbon. Thus, the effective 
gap of the junction decreases as the length increases. As 
expected, the transmission function of the pristine infi- 
nite 3-aGNR acts like an envelope curve for the other 
curves; it has an energy gap of ~1.5 eVj^ This aspect 
confirms that the contacts between the graphene sheets 
and the nanoribbon do not represent significant barriers 
to the transport of electrons, as one can naively expect 
due to the chemical identity of the different subsystems. 
Within the energy gap of the nanoribbon, the transmis- 
sion function is very low and decreases rapidly as the 
length of the ribbon increases. This behavior is clearly 
due to the tunneling mechanism dominating the trans- 
mission at those energies. The states of the graphene 
lead decay exponentially along the ribbon because there 
are no states in the junction supporting the conductance. 
This metal-semiconductor-metal device presents an effec- 
tive gap for transport, since the conductance inside the 
gap is several orders of magnitude lower than outside. 

We now discuss the effects of the aGNR width in the 
transmission function, considering junctions made with 
4L 3-aGNR, 4L 5-aGNR, and 4L 7-aGNR. We consider 
a representative aGNR for each of the three classes, so 
that we can capture the main differences among them. 
As illustrated in [SJ the aGNR junctions have different 
gaps in the transmission spectrum consistently with the 



previously reported results J^ In miniaturized graphene- 
based electronic devices these structures may act both as 
linkers and as active components, and it would be use- 
ful to finely tune their width in order to have different 
conductive behaviors. Within each class, we expect the 
energy gap to decrease as the width of the junction in- 
creases. We notice that the intensity of the transmission 
function is higher for wider ribbons, as more electronic 
channels are open. Moreover, as the width of the rib- 
bon increases the spike features of the transmission func- 
tion become less marked, since the system approaches the 
limit of graphene. In this case, the transmission functions 
of the pristine aGNR is a good reference only for those 
junctions where the aspect ratio between the length and 
the width is reasonably high. This because the contri- 
bution to the conductance of the infinite aGNR comes 
from fc|| = only. For finite nanoribbons this restriction 
is relaxed and the passage of electrons with a finite fcii 
is allowed, with fcii within an interval that increases as 
the width to length ratio increases. This is clearly re- 
flected in the transmission function, which departs from 
an exponential decay within the energy gap for short and 
wide ribbons. This is more evident for the 4L 7-aGNR 
in [5l For this ribbon the conductance is linear in the 
gap region, reflecting the underlying electronic structure 
of the graphene electrodes, and departing from a strictly 
tunneling behavior dominated by the energy dependence 
of the effective tunneling barrier. 

We now discuss the effects of the application of a bias 
to the junction. We consider the junction with the 8L 
3-aGNR as a representative case; this is the longest junc- 
tion considered here, thus being that in which the edges 
of the two graphene leads are less interacting. In [5] the 
transmission functions for different biases up to 1.0 V are 
shown. For any given bias voltage, the chemical poten- 
tials of the two leads are well recognizable, as they corre- 
spond to the left and right Dirac points where the DOS 
tends to zero. Within the bias window, the intensity of 
the transmission function increases with the bias. This 
is expected since it is determined by the number of elec- 
tron tunneling from the leads, and the DOS of graphene 
shows a linear energy dependence in that region. Thus, 
the shape of the transmission function within the gap re- 
gion can be interpreted in terms of the product of the 
DOS of the left an right graphene leads and a modula- 
tion function determined by the size of the gap of the 
nanoribbons. We notice two small peaks near the Dirac 
points of the two electrodes, which we assign to the edge 
states appearing at the two zigzag electrodes42i In fact, in 
presence of zigzag edges there is a band near Ep which 
corresponds to a state localized at the edge. At small 
applied bias, the probability for an electron lying in an 
edge state to be transferred to the opposite lead is very 
small since the DOS at that energy is negligible. How- 
ever, at finite bias, the situation is different and we can 
find distinct features in the transmission curves related 
to the presence of edge states. We can also sec that the 
resonances originated by the HOMO and LUMO of the 
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FIG. 6. Lower panel: current-voltage characteristics of 8L 
3-aGNR junction. Upper panels: transmission function of 8L 
3-aGNR junction for four selected values of the bias voltage, 
as indicated in labels. The bias window is delimited by two 
vertical lines. 



aGNR slightly shift as the bias is applied. As a result, 
the energy gap reduces by approximatively 10% when a 
bias of 1.0 V is applied. We verified that the current re- 
mains very low (less that 0.1 ^A) for biases up to 1.0 V. 
For larger values, the intensity sharply increases when 
the two main peaks enter the bias window, and for a bias 
of 3.0 V we calculated a current of 3.3 ^A. 

When considering junctions made of 3-aGNR, the sys- 
tem is actually a chain of phenyl rings. In this case, it 
has been shown that the stable configuration consists of 
neighbouring rings lying on different planes^i*^ We can 
thus take into account one more degree of freedom, that 
is the relative torsion angle between the phenyl rings of 
the linkers. We show the effect of the ring rotation by 
taking the 4L 3-aGNR as test system. According to our 
calculations, the most stable configuration corresponds 
to a rotation of ~ 45°4^ In fact, the central ring tends 
to rotate with respect to its planar configuration due to 
steric repulsion between its hydrogen atoms and those of 




FIG. 7. Transmission function at zero bias for the 4L 3-aGNR 
junction, under rigid torsion of the central phenyl ring. The 
curves shown here refer to rotations by 0°, 30°, 60°, and 90°. 
The main resonances are highlighted and the labels corre- 
spond to the 3-aGNR molecular states shown in [51 



the neighboring rings. By looking at the transmission 
function in [7] we can note a general trend. As the ring 
rotates from 0° to 90°, the peak H (HOMO) shifts to- 
wards lower energies, while the peak L (LUMO) shifts 
towards higher energies. As a result, the energy gap 
of the system increases. In fact, the wavefunctions re- 
lated to these peaks (see [8|) are extended along all the 
molecular backbone and they feel a distortion as the cen- 
tral ring is rotated, thus their energies are expected to 
change. Both HOMO and LUMO are even respect to 
the mirror plane that bisects the 3-aGNR molecule. As 
consequence of this symmetry, when the rotation is by 
90° these two peaks disappear: the linear combinations 
of the Pz orbitals of the central carbon atoms are orthog- 
onal to those of the remaining carbon atoms, resulting 
in the closure of the corresponding conduction channels. 
The situation is different for HOMO-1 and LUMO-f 1, 
whose wavefunctions are mainly localized on the central 
phenyl ring. In fact, they are at the same energy position 
regardless of the torsion angle. These considerations may 
be extended to the case of longer aGNR. 

Concerning the passivation of the graphene edges, all 
the results presented so far refer to what we can call 
"standard saturation". That is, each edge carbon atom 
is saturated by a single hydrogen (panel a in [9]). In a real 
system subject to different hydrogen partial pressures, 
one may guess that the edge saturation can change4S In 
order to validate our results, we show the effect of four 
different edge saturations on the transmission function. 
The first case we consider is the formation of pentagon- 
heptagon reconstruction of the graphene edges, where no 
H atoms are present on the edge (panel d). Then, we 
considered the case in which each carbon edge atom is 
saturated by two hydrogens (panel b). One may also 
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chemical species constituting the subsystems. Hence the 
transport properties are mainly determined by the shape 
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FIG. 8. Isosurface plots of the selected molecular orbitals 
of the isolated 3-aGNR linker: HOMO (H), HOMO-1 (H-1), 
LUMO (L) and LUMO+1 (L+1). 



guess a more complicated connection between the aGNR 
and the leads, and we modeled it as shown in the inset of 
panel c. [9] shows the transmission function for these four 
cases. Our results are very robust with respect to these 
changes of the geometry and saturation of the graphcne 
lead edges, as we see only minor differences among these 
cases, mainly a sharpening of the peaks derived from the 
HOMO and LUMO resonances. This is evident in par- 
ticular for the double-H and shaped edges, in which the 
orbitals of the aGNR are less hybridized with the leads. 
The reason is that in both cases less charge is present on 
the graphcne layer edge. 



IV. CONCLUSIONS 

Very simple prototypes of carbon-based molecular 
junctions have been investigated in order to character- 
ize their transport properties. Carbon nanoribbons with 
armchair shaped edges have been considered as linkers 
between between two semiinfinite graphene electrodes. 
These systems present a typical metal-semiconductor- 
metal behavior due to the electronic gap of the ribbons 
which depends on the their width. The electronic prop- 
erties of the isolated subsystems, i.e. graphene electrodes 
and the nanoribbon, together with the interaction in con- 
tact regions determine the transport properties of the 
junction. For what concerns the coupling between the 
subsystems, this is very efficient and the contacts do not 
create appreciable barrier for transport due to the same 







FIG. 9. Zero bias transmission function of 8L 3-aGNR junc- 
tions with four different saturations of the graphene edges. 
The inset show the linking geometries. 



(width and length) of the finite ribbon included in be- 
tween the graphcne leads. Larger molecules furnish more 
channels for convoying the electrons then the smaller ones 
and consequently the conductance is generally higher. In 
the energy gap region the transport occurs via electron 
tunneling between the electrodes and the efficiency de- 
cays simply with the linker length. 

Some additional representative configurations has been 
considered explicitly: thinest ribbon a torqucd geome- 
tries has been studied to evaluate the effects of the mis- 
alignment of the TT orbitals on different phenyl rings. The 
changes in transmission through HOMO and LUMO or- 
bitals are noteworthy, while for energy levels farther from 
Ep they are negligible. Finally, graphene edges with dif- 
ferent saturations and reconstructions has been consid- 
ered and they do not influence significantly the calcu- 
lated transport properties. In spite of the simple model 
adopted in this paper, the reported results on the elec- 
tronic transport provide some parameters to control and 
engineer future carbon-based electronic devices. 
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